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Introduction Cold atomic gases in optical 
lattices are emerging as excellent laborato- 
ries for testing models of strongly interacting 
particles in condensed matter physics. Cur- 
rently, one of the major open questions is how 
to obtain the finite temperature phase dia- 
gram of a given quantum Hamiltonian directly 
from experiments. Previous work in this di- 
rection required quantum Monte Carlo sim- 
ulations to directly model the experimental 
situation in order to extract quantitative in- 
formation, clearly defeating the purpose of an 
optical lattice emulator. Here we propose a 
new method that utilizes deviations from a 
local fluctuation dissipation theorem to con- 
struct a finite temperature phase diagram, for 
the first time, from local observables accessi- 
ble by in situ experimental observations. Our 
approach extends the utility of the fluctuation- 
dissipation theorem from thermometry to the 
identification of quantum phases, associated 
energy scales and the quantum critical region. 
We test our ideas using state-of-the-art large- 
scale quantum Monte Carlo simulations of the 
two-dimensional Bose Hubbard model. 

As thermal fluctuations decrease, interactions be- 
tween particles dominate and drive systems into 
emergent quantum phases of matter such as super- 
fluids, Mott insulators, magnetically ordered phases, 
and spin liquids [1, 2, 3, 4]. Not only are these phases 
with their associated excitations interesting, but the 
transition between quantum phases opens up a quan- 
tum critical region dominated by large fluctuations. 
The large fluctuations arise from new degrees of free- 
dom that must form as the system transits from one 
phase to the other, tuned by a coupling parameter as 
illustrated in Fig. 1. 



Quantum gases confined in an optical lattice of- 
fer a unique platform to study such quantum phase 
transitions since the coupling strength is easily tuned 
by the laser intensity or the magnetic field. Detec- 
tion and characterization of these systems has, until 
recently, been limited to time-of-flight observations. 
Experiments on bosons in optical lattices, emulating 
the Bose Hubbard model (BHM), have shown signa- 
tures of the coherent superfluid (sharp peaks) and 
of the Mott state (broad, featureless distributions) 
[1, 5]. Further theoretical investigation has indicated 
that sharp peaks in the momentum distribution can 
occur even in the non-superfluid phase and there is a 
wealth of information about quantum critical fluctu- 
ations in the nature of these peaks [6, 7]. 

The primary bottleneck in obtaining a more quan- 
titative comparison with theory has proven to be a 
definitive measurement of the temperature and a di- 
agnostic of the all-important quantum critical regime. 
Some preliminary identifications have been made, but 
they require input from individually tailored large- 
scale quantum Monte Carlo (QMC) simulations [8, 9]. 

In this Letter, we provide the theoretical frame- 
work for the application of a local fluctuation dissi- 
pation (LFD) theorem to experiments to gain funda- 
mental insights into the nature of the phases, their 
low-lying excitations and quantum criticality. Moti- 
vated by the success of recent experiments [10, 11, 12, 
13, 14] to access the local density in situ, we define a 
LFD ratio R as the ratio of the local compressibility 
Ki to local number fluctuations 5nf at site i, defined 

by 

R = Ki/5n\ (1) 
n t = d{ni)/dn 

Remarkably, from the temperature dependence of 
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Figure 1: Identifying Phases and the Quantum Critical Region using Deviations from the Fluctuation Dissi- 
pation Theorem: (a) Finite temperature phase diagram of a quantum Hamiltonian, e.g. the Bose Hubbard 
model, shows the suppression of the superfiuid transition temperature T c upon repulsive interaction U in- 
creasing relative to tunneling t and ultimately vanishing at a critical coupling (t/U) c . For t/U < (t/U) c , 
the system is Mott insulating with a crossover scale Xa that also vanishes at (t/U) c . T* is the temperature 
scale at which local quantum effects begin to emerge. The shaded region above (t/U) c depicts the quantum 
critical region dominated by large fluctuations, (b) We show how such a phase diagram for a given quantum 
Hamiltonian can be obtained from the local fluctuation-dissipation ratio R defined in equation (2). In the 
gapped Mott phase we extract /3 ma x, which is related to the inverse Ta as discussed in the text; similarly, 
the maximum in the gapless SF phase /3 max corresponds to the inverse T c . For both phases the deviation of 
R from linearity marks (3* = 1/T* where quantum effects begin. In the critical region, R remains linear in 
/3 but with a proportionality constant less than unity. 



this single quantity R, we propose that it is possible 
to estimate the temperature of the onset of quan- 
tum effects T* , the temperature of quantum phase 
ordering T max and the quantum critical region, for a 
general quantum Hamiltonian. We test our proposal 
on the 2D BHM using large-scale QMC simulations 
and summarize our central results for phase identi- 
fication in Fig. 1(b) and in the finite-temperature 
phase diagram in Fig. 4, obtained solely from this 
LFD ratio. This explicit demonstration opens up the 
possibility of finding such phase diagrams for gen- 
eral quantum Hamiltonians directly from experimen- 
tal data. We thus provide the crucial missing link 
in the grand challenge to emulate strongly correlated 
materials such as the high temperature superconduc- 
tors using ultracold atoms in optical lattices. 
Fluctuation-Dissipation Theorem (FDT): The 
full quantum FDT relates the imaginary response 
function x" to the dynamic structure factor S at in- 
verse temperature /3 and is given by 
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X"(q > w)= S(q,w) 



(2) 



Here we specialize to perturbations of the density 

;k a k+q C 



Pil) = Sk a ik+a a k arL d the density-density corre- 



lation function 5(q, t) = (p(q, t)p(0, 0)). A con- 
served quantity, such as the total number of parti- 
cles N = p(q = 0), commutes with the Hamiltonian, 
so p(q = 0, t) is independent of t. As discussed in 
the supplement, this leads to the exact expression 
X(q -> 0,cj = 0) = (5 S(q -> 0,t = 0). The zero 
frequency response function equals dnjd^i = t?kt 
where kt is the thermodynamic compressibility, and 
the equal-time correlation function equals (3 5N 2 /V 
where SN 2 = (N 2 ) — (N) 2 is the total number fluc- 
tuations. Equating the two gives 



2 

n kt 



d(n) _ 5N 2 



Of i 



V 



(3) 



While equation (3) resembles the classical FDT, it is 
crucial to note that it is valid even in the quantum 
regime since p(q = 0) is a conserved quantity. 

There are two ways to probe the system locally. We 
could locally perturb the system with a small change 
in the chemical potential 5/ii at site i and measure 
the resulting local density variation Srii. This pro- 
cedure leads to the single site form of the classical 
FDT, Ki l = = (3 Snf, that is only valid at high 
temperatures since the density on a given site is not 
conserved. Both forms of the FDT relations in equa- 
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tion (3) and in the single site form Ki l discussed 
above enable thermometry [15, 16, 9], but serve a 
limited purpose in revealing other properties of the 
system. 

We propose that a more useful quantity, sensitive 
to the nature of the phase and to various energy 
scales, is the local fluctuation-dissipation (LFD) ratio 
R defined in equation (2). The essential difference is 
that R involves Kj = d(rii) /dfi which is the measured 
change of the local density in response to a global 
chemical potential [i variation, and as such is sensitive 
to long range order and phase transitions, as opposed 
to Ki l which is the local response to a local pertur- 
bation. For a homogeneous system, «j is the same as 
the thermodynamic compressibility (up to factors of 
n 2 ). Although we investigate R in a uniform system 
below, it is important to keep in mind the applicabil- 
ity to inhomogeneous systems as created in trapped 
atomic gases. Assuming local density approximation, 
Ki and R can be extracted directly from the density 
profile n(n(r)) [17]. 

In the following, we demonstrate the potential for 
R to map out finite temperature phase diagrams by 
evaluating R in the BHM. Bosons trapped in an op- 
tical lattice and confined in a potential are described 
by 

H = - - ^2(a\aj + Oiot) 

Z (i,3) 

+ ~'^2n i (n i -l)-'^2nin i . (4) 

i i 

Here dj (<z|) are boson annihilation (creation) op- 
erators, (. . .) indicates nearest neighbor sites, \ii = 
/io — ctrf is the chemical potential on site i for a 
parabolic confining potential, and t and U set the 
hopping and interaction energy scales, respectively. 
We simulate the BHM at finite temperatures with 
worldline QMC using the directed loop algorithm on 
up to 32 2 site lattices [18]. We establish the essential 
ideas in a uniform system (a = 0) for clarity, but the 
results are easily extended to the nonuniform system 
by using the local density approximation. 
Classical Regime: At high temperatures, T 3> 
U ^> t, the system is in a non-interacting classical 
regime. Here the inverse LFD ratio R^ 1 sa T is de- 
termined primarily by the density fluctuations dnf 
that increase linearly with temperature («j remains 
finite and independent of T). As T decreases to the 
regime U 3> T 3> t, the system remains classical. 
Interaction effects cause both Ki and Snf to deviate 
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Figure 2: High Temperature Behavior of Local 
Observables. (a) The local compressibility Ki = 
d{ni)/d^ji and the local number fluctuations Snf well 
inside the SF phase of the 2D BHM. As expected for 
an ideal gas, k, approaches a constant and Snf grows 
with T at high temperatures, (b) R^ 1 = 5nf/rt, ~ T 
confirms the high temperature behavior (dashed line) 
of the LFD ratio, equation (2). Deviation of J? -1 
from linear behavior allows us to determine the onset 
of quantum effects. 

from the high-T ideal gas limit but the inverse LFD 
ratio remains approximately linear in T (see Fig. 2). 

R can also be used to test for equilibration of the 
system in different regions. For bosons in optical lat- 
tices, the Mott-like center is sometimes observed to 
be at a significantly different temperature from the 
superfluid or normal wings [19]. In such situations, 
successive local measurements of R can be useful to 
garner information about rate-limiting processes for 
achieving equilibrium. 

Onset of Quantum Effects T*: The deviation of 
R^ 1 from linear T behavior defines the temperature 
T* at which quantum effects first become evident. 
For bosonic systems at low density n < 0.1 we observe 
bunching (see supplement Fig. 1), an enhancement of 
occupancy at a site, due to quantum statistics and 
manifested as R> (3. 
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Figure 3: Phase Signatures and Energy Scales Encoded in R. Panel (e) The T = BHM phase diagram 
with the parameters of each dataset are indicated. Panels (a) and (b): In the MI, the compressibility Ki 
(circles) is suppressed at low temperatures to zero while local quantum fluctuations maintain Snf (squares) 
at a finite value. The corresponding LFD R — Ki/5nf in the MI shows a characteristic maximum and 
decays exponentially at low T, from which we determine the particle-hole gap energy A p h discussed in the 
Supplement. Panels (c) and (d): In the SF, Ki, Snf and R all approach a finite constant as T — > 0. The peak 
in R occurs at f3 max and agrees well with the inverse superfluid T c (indicated by the arrow) obtained from 
the superfluid density p s (squares in (d)). In all panels (a,b,c,d), the inverse temperature (i* indicates when 
R deviates significantly from the classical limit. The dashed black line indicates the high temperature R sa (3 
limit in (b,d) and horizontal bars indicate uncertainty in f3* and f3 m ax in each panel. Panel (f) shows the 
behavior close to the quantum critical point, where the system does not order at the temperatures probed. 
Both Snf and Ki (not shown) do not converge to the behavior of either phase, and R oc /3 with a slope less 
than one indicates the absence of an energy scale, a characteristic of the quantum critical region. 



At higher density, the bunching tendency 
from statistics competes with inter-boson re- 
pulsion that tends to keep bosons apart. In 
this regime, density correlations between sites 
are anti-correlated (see supplement Fig. 2). It 
can be shown in general that R = Ki/Snf — 

[Ei^ - («.•) /(VSnf)) which 

implies that for anti-correlated density fluctuations, 
R < (3 as seen in Fig. 2. 

Phases: Understanding the behavior of Ki and Snf, 
the observables composing the LFD ratio, is essential 
for understanding the signatures of the phases in R 
(see Fig. 3 for typical MI and SF systems). The peak 
in R at a temperature T max is a generic feature of the 
system entering an ordered phase. 

Superfluid T c : The gapless collective excitations in 
the SF cause both Ki and Snf to approach a constant 



value as T — ► 0. As the system condenses at T c , the 
compressibility Ki, a long wavelength response func- 
tion, exhibits a peak because of critical fluctuations. 
Since the appearance of long-range phase coherence 
does not affect the smoothly decreasing local Snf, the 
peak in R mirrors the peak in Ki near T c . Compari- 
son with the SF T c calculated directly from our QMC 
simulations, determined by the vanishing of the su- 
perfluid density, confirms that T max « T c , illustrated 
in Fig. 3(d). 

Mott Insulator Ta: The Mott gap suppresses the 
low energy excitations contributing to Ki, causing Ki 
to vanish as T — > 0. On the other hand, in spite 
of the Mott gap, the local number fluctuations re- 
main finite down to the lowest temperatures because 
of local quantum fluctuations. As the temperature is 
increased, in contrast to the SF-normal phase transi- 
tion, the MI crosses over into the normal state with no 
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Figure 4: Determination of the Phase Diagram of a Quantum Hamiltonian using Local Observables: This 
figure shows the contour plot of LFD ratio R in the inverse temperature-chemical potential {fit — li/U) 
plane at (a) t/U — 0.15 and (b) t/U — 0.4. We determine the phase diagram from R{[3) by marking 
two inverse temperature scales, (3* = 1/T* (solid squares) that indicates the onset of quantum effects, 
and fSmax = 1/T max (solid circles) that delineate phase boundaries or crossovers. (3 m ax agrees with the 
characteristic temperature scales in the superfluid and the MI phases, or 1/T C (open circles) and 1/Ta (open 
squares), respectively. Here T c and Ta are calculated independently from R by using the vanishing of p s 
and ft, and confirm the interpretation of the T max scale. As seen from Fig. 3(a), at T = and t/U = 0.15, 
the system progresses from vacuum to SF, hits two quantum critical points bounding the MI, and returns 
to SF as li/U increases in (a), while in (b) it moves from vacuum to SF and crosses no critical points. The 
signatures in R of each state, identified in Figs. 2 and 3, are visible here as a linear increase in R at high 
T, a plateau in the SF and a decay to zero in the MI. In panel (b), note that even with no MI, the window 
of critical fluctuations at temperatures T > T c of the strongly coupled SF varies with /i/U because of the 
proximity of the t/U ~ 0.24 quantum critical point. 



transition and therefore shows no specific signature in 
k. Thus, the peaks of R and Ki do not necessarily line 
up. We identify the peak in R with the MI crossover 
temperature Ta as confirmed again by QMC simu- 
lations, where Ta corresponds to vanishingly small 
compressibility (see Fig. 4(a,b)). We further deter- 
mine the zero temperature energy gap to add 
a particle or a hole, whichever is smaller [20] by fit- 
ting Ki or R by e _/3Aph within the MI. The extracted 
A p h agrees very well with QMC simulations [21] deep 
in the MI (Fig. 3(e)), but differ on approaching the 
critical (fi/U) c points, which we attribute mainly to 
finite temperature and size effects. 

We put together our knowledge of the signatures 
in R in the various quantum and classical phases to 
construct the phase diagram in Figs. 4(a,b) at cou- 
plings above and below (t/U) c . The density in each 
diagram changes from vacuum to n ~ 1.5 particles 
per site. 



Critical Regime: The degeneracy temperature T* 
depends on t/U and fi/U and is reassuringly indepen- 
dent of the underlying critical points (see Fig. 4(a)). 
In the temperature range between T* and T max , 
quantum critical fluctuations lead to R cx (3 with a 
slope distinctly less than unity (Fig. 4(e)), and Ki and 
Snf also display nontrivial behavior. This region is 
clearly largest near the critical points in Fig. 4(a), but 
the proximity of the critical coupling (t/U) c rs 0.24 
opens up a quantum critical region at intermediate 
density and temperature in Fig. 4(b) as well. For 
both t/U couplings, the window of critical behavior 
narrows as the density increases, reflecting a change 
from quantum criticality (proximity to quantum crit- 
ical point) to classical criticality (as expected near 
the SF-normal phase transition). 

In conclusion, while both the thermodynamic 
dN/dLL = (3 SN 2 and the local dru/dfii = /3 5n 2 FDT 
are exact relations useful for estimating the tempera- 
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ture, what is new in our proposal is the construction 
of a LFD ratio R defined in equation (2), that involves 
measuring local density fluctuations in response to a 
global chemical potential change. We show that R 
is sensitive to far more than just the temperature. 
As proof of principle, we have demonstrated that R 
identifies phases and critical regimes as well as esti- 
mates the Mott mobility gap A p /j. It should therefore 
be possible to experimentally map out finite temper- 
ature phase diagrams, as in Fig. 4(a,b), without the 
need for individualized QMC simulations. It is also 
possible to extend R to other quantities such as prob- 
ing spin susceptibility and corresponding spin fluctu- 
ations for magnetic systems. Given the very funda- 
mental basis on which the LFD ratio is constructed, 
we expect it to be an ideal candidate for probing 
phases and quantum criticality of general quantum 
Hamiltonians. 
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1 Bunching or Anti-Bunching 

The quantity R depends on the ratio of the nonlocal number fluctuations, SfiiSrij = (niUj) — 
(rii) (rij) between sites i and j, to the local number fluctuations 5nf = 8ni5rii. Specifically, 

At high temperatures, the nonlocal number fluctuations are vanishingly small and R ~ /3. 
Any finite nonlocal number fluctuations cause R to deviate from /3. Boson statistics favor 
bunching, or positively correlated nonlocal number fluctuations, and this bunching behavior 
is observed in R in the Bose-Hubbard model at very low densities. Bunching implies the 
second term in equation (1) is positive and R > /3, as shown in Fig. 1. 




2/3* 6 

Figure 1: Boson Bunching due to Statistics at Low Density. The main plot shows R (circles) 
becomes larger than (dashed line) at low temperatures. The inset plots the density (circles) 
over the same temperature range. The parameters are t/U = 0.15, fi/U = —0.2 and the error 
in 0* is indicated by the horizontal bar. Note that this system does not become superfluid 
at the temperatures probed. 

Anti-bunching arises from the repulsive interactions in the BHM. Interactions energeti- 
cally favor negatively correlated nonlocal number fluctuations and overwhelm the statistical 



tendency to bunch at finite density (> 0.1) and at (t/U) near the critical (t/U) c . This 
implies the second term in equation (1) is negative and R < (3 at low temperatures. This is 
observed to be the case in any quantum phase near (t/U) c and is exemplified in Fig. 2. 




Figure 2: Interactions cause Bosons to Anti-Bunch at Moderate Densities. The main plot 
shows R (circles) becomes smaller than (5 (dashed line) at low temperatures when the system 
enters a quantum phase (superfluid in this case). The large error in f3 max (horizontal bar) is 
characteristic of the superfluid in the presence of large interactions, where the peak in R is 
difficult to resolve. The inset shows density (circles) as a function of /3. The parameters are 
t/U = 0.15, n/U = 0.05. 
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2 The Particle-Hole Gap 



The particle-hole gap A p h is the minimal energy required to either insert or remove a particle 
from a Mott state. This energy scale can be extracted from R by fitting its decay from the 
peak to zero with the simple exponential form e~^ Aph . This form also describes the decay of 
the compressibility Ki as a function of /3 within the Mott state. We compare extracted A p h 
with T = QMC results from Capogrosso-Sansone, et al, in Fig. 3. 
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Figure 3: Extracting the Particle-Hole Gap from R. The righthand column shows R (circles) 
at various mu/U values and the exponential fit to extract A p h (solid line). Plotted in the 
lefthand figure is the extracted A p h (circles) and the results from T = QMC calculations 
by Capogrosso-Sansone, et al. 



3 The Fluctuation-Dissipation Theorem 

The fluctuation-dissipation theorem (FDT) relates the imaginary part of the generalized sus- 
ceptibility to the correlation function. For completeness, we discuss (a) the salient points in 
the derivation of the quantum FDT; (b) the static structure factor; (c) the high temperature 
limit of the quantum FDT and its reduction to the classical FDT. (d) Further, we show that 
for a conserved quantity, the classical FDT expression is valid even in the quantum regime, 
which may seem somewhat surprising. And lastly, as a test, we provide the compressibility 
and the number fluctuations of the ideal gas in the high temperature limit. 

3.1 The Quantum Fluctuation-Dissipation Theorem 

Consider a quantum system of volume Q defined by a time-independent Hamiltonian H with 
many body states and energies H |^ n ) = e n \^Sf n ), perturbed by a probe H'(t). We assume 
that the external perturbation F^t) couples to an operator A of the system via H'(t) = 
-AF A {t). For spatially varying probe H'\t) = - J drA(r)F A (r,t) = £ q A_ q F q (t). For 
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example, in our case A(r) is the density operator n(r) = ^(^^(r) = ^ e iq ' r 7i(q) where 

™(q) = Ek«Lq/2 a k+q/2. 

The response (B) to linear order in the perturbation is (B)(t) = dt'xBA(t — t')FA(f) 

where Xba{P — t') = iO(t — t') ^ B(t), A(t')J ^. By using a spectral representation in terms of 

exact eigenstates of H and the Heisenberg representation of the time dependent operators, 
we obtain 

(^— q)mn(-^q)nm (■^q)mn (^— q)nr 



1 v-^ e - ^ 



w + «?7 + e„ m u + irj- e r 



(2) 



where e ram = e n — e m and rj = + is a small positive number to ensure convergence as 
t — > oo. The well-known identity lim = P(— ) ± 27r<5(:r) yields the imaginary part of 

»7^0+ X ±11] X 

the response function, 

XsA(q> W ) = q ^ [(^-q)mn(-Bq)nm5(o; - € nm ) - (B^ mn (A_^) nm 8(u + 6 nm )] . (3) 

m,n 

Next we consider the corresponding correlation function typically measured in a scattering 
experiment defined by 

S BA (r,t;r',t') = (B(r,t)A(r',t')) (4) 

5 flA (q,t-0 = ^(B q (t)i- q (0) (5) 
for a translationally invariant system. Using the spectral representation, we obtain 

27T . g-Z^m 

^(q, w) = -jy — ^— (5 q ) m „(A_ q )„ m 5(a; - e nm ). (6) 

m,n 

By exchanging the indices in the second term in Eq.3, we obtain the quantum fluctuation 
dissipation theorem (QFDT) 

XBA(<liU) = T^ 1 " e "^) -^( B q)"»(^-q)™^(w - "urn) (7) 



1 - e-P" 



2 

3.1.1 Static Structure Factor 

The static structure factor is defined by 



Sba(cl,ou). (8) 



Sba(<1) = S BA (q,t = 0) = / ^^(q,^) (9) 

= Lsi=FS^<* w) - < 10 > 
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Using the oddness property x'ba(~ u ) = ~Xba( uj ) yields 

Sba(<i) = ± (B q i_ q ) = /°° v coth (^r) ( n ) 

3.1.2 The High Temperature Limit of the QFDT 

At temperatures k B T ^> hu larger than any characteristic frequencies of the system, coth (4jr) — > 
2/(3u and the static structure factor reduces to 

P Jo 71 CO 

= k B T X ' BA (c l ,cu = 0) (13) 



where we have used the Kramers-Kronig relation x' BA (q, u) = P ~ Xb ^ ~ to relate 
the real and imaginary parts of the response function. Since at cu = the imaginary part is 
zero, we can replace x' by simply x- 



3.2 The QFDT for Conserved Quantities 



We next use (i) the definition of the correlation function for a conserved quantity, (ii) the 
quantum FDT, and (iii) Kramers-Kronig relation to finally derive Xaa(q — > 0,w = 0) = 
P~Saa{<I — > 0,t — 0). The derivation is detailed below. 

A conserved quantity A(q = 0) = A such as n(q = 0) = ^ k a^a^ = N, the total number 



of particles, commutes with the Hamiltonian, 



An,H 



0. This implies that the matrix 



element (\I/ r 



A ,H 



ty n ) = or equivalently (e n — e m )(Ao) m , n = 0. If m ^ n, we must 



have (A ) mtn = which results in 



lim lim Xaa((1, u) = lim — - 

m,n 



-Pen 



(A)) 



12 1 



(u + it]) 2 - e 



2 

nrn 







(14) 



using Eq. 2. Thus Xaa(q = 0, u — > 0) = if A(q_ = 0) is a conserved quantity. 
However, the situation is totally different if we change the order of limits. 



lim lim Xaa(<1, w) = lim — > 



-Per, 



^nrn | (^q)mn 

(u + it]) 2 - e 2 



n 2 K 



(15) 



For density operators, the last equality in Eq. 15 follows from the perturbation H' = 
— J dr8n(r)8fi(r,t) = — ^ ^^Sfi-qSji^t) which produces a response in the system 



Xnn(q -> 0,W = 0) 



dn 



-(to = 0) 



n 2 K 



(16) 
(17) 
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where k is the isothermal compressibility. 

Another way to understand the behavior of QFDT for a conserved quantity (^A(q = 0, t) 



is to note that it is independent of t. In addition the correlator \^A(q = 0,t)A(— q = 0, t') 

is independent oit — t' and hence its Fourier transform must be a delta function in frequency. 
Thus from Eq. 5 

S AA (q = 0, u) = 2nS(uj)S AA (q = 0) (18) 

From the quantum FDT Eq. 8 we obtain for a conserved quantity whose correlation 
function is given by Eq. 18, 

x'Uq = 0, w) = (1 - e-^)nS AA (q = 0)6 (u) (19) 
By using the Kramers- Kronig relation we get 

XAA (q = 0,u^0) = S AA (q = 0) / duo 8{u) (20) 

= (lS AA (q = 0) = !L (P> 

For density operators we thus get 



(5S AA (q = Q)='-^A\q = Q)) (21) 



Xnn(q = 0) = ^<n 2 (q = 0))=n 2 K (22) 

We would like to stress that while this result may look similar to the high temperature limit 
of equation (13), it is valid in the quantum regime for a conserved quantity. 

3.3 The Ideal Gas 

We derive the high temperature behavior of the compressibility, dn/dfi, and number fluctu- 
ations 5N 2 = (N 2 ) — (N) 2 of the ideal gas. From the equation of state PV = Nk B T, we 
get 

1 (dV\ 1 On _ (3 

KT = ~V\dP)^ N = ^-n (23) 

where n = N/V. 

For an ideal gas, the chemical potential is related to the density by [3 fx = — log ^^rj, 

where the thermal deBroglie wavelength At = h/ \f2nmk,BT . For a fixed fx, the high tem- 
perature expansion of n(/i, T) is 



which implies that the temperature dependence of the local compressibility is 

n 2 «=^~T d / 2 - 1 (l + Ai T- 1 ). (25) 

OIL \ / 



2" 
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Using the local fluctuation-dissipation theorem at high temperatures, we find 

( < Jn 2 >»^r~r d / 2 (l + A xT- 1 ), (26) 

as we might have guessed from a simple Brownian motion or diffusion model of the number 
fluctuations. 

Thus, we find that in 2D, while the local compressibility n 2 n is independent of tempera- 
ture (note that in the Letter we have absorbed the n 2 factor in the definition of k) and the 
number fluctuations (5n 2 ) scale linearly with T . 

In the presence of interactions, both number fluctuations and compressibility deviate 
from their classical values, however their ratio R = u 2 k/ (5n 2 ) continues to be linear in 
T _1 as determined by the fluctuation-dissiaption theorem (see Fig. 2 of Letter). Only when 
quantum effects become important does R exhibit non-trivial behavior. 



7 



